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WASP-17b: AN ULTRA-LOW DENSITY PLANET IN A PROBABLE RETROGRADE ORBITt 
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ABSTRACT 

We report the discovery of the transiting giant planet WASP-17b, the least-dense planet currently 
known. It is 1.6 Saturn masses but 1.5-2 Jupiter radii, giving a density of 6-14 per cent that of Jupiter. 
WASP-17b is in a 3.7-day orbit around a sub-solar metallicity, V = 11.6, F6 star. Preliminary detection 
of the Rossiter-McLaughlin effect suggests that WASP-17b is in a retrograde orbit (A « — ISOdeg), 
indicative of a violent history involving planet-planet or star-planet scattering. 

WASP-17b's bloated radius could be due to tidal heating resulting from recent or ongoing tidal 
circularisation of an eccentric orbit, such as the highly eccentric orbits that typically result from 
scattering interactions. It will thus be important to determine more precisely the current orbital 
eccentricity by further high-precision radial velocity measurements or by timing the secondary eclipse, 
both to reduce the uncertainty on the planet's radius and to test tidal-heating models. Owing to its 
low surface gravity, WASP-17b's atmosphere has the largest scale height of any known planet, making 
it a good target for transmission spectroscopy. 

Subject headings: planetary systems: individual: WASP-17b - stars: individual: WASP-17 



1. INTRODUCTION 

The first measurement of the radius and density of 
an extrasolar planet was made when HD 209458b was 
seen to transit its parent star (Charbonneau et al. 2000, 
Henry et al. 2000). The large radius (1.32 Rjup) of 
HD 209458b, confirmed by later observations (e.g., Knut- 
son et al. 2007), could not be explained by standard 
models of planet evolution (Guillot & Showman 2002). 
Since the discovery of HD 209458b, other bloated plan- 
ets have been found, including TrES-4 (Mandushev et al. 
2007), WASP-12b (Hebb et al. 2008), WASP-4b (Wil- 
son et al. 2008; Gillon et al. 2009a, Winn et al. 2009a, 
Southworth et al. 2009), WASP-6b (Gillon et al. 2009b), 
X0-3b (Johns-KruU et al. 2008; Winn et al. 2008) and 
HAT-P-lb (Bakos et al. 2007; Winn et al. 2007; Johnson 
et al. 2008). Of those, TrES-4 is the most bloated, with 
a density 15 per cent that of Jupiter, and a radius larger 
by a factor 1.78 (Sozzetti et al. 2009). 

The mass, composition and evolution history of a 
planet determines its current radius (e.g.. Burrows et al. 
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2007; Fortney et al. 2007). Recently, numerous theo- 
retical studies have attempted to discover the reasons 
why some short-orbit, giant planets are bloated. A small 
fraction of stellar insolation energy would be sufficient to 
account for bloating, but no known mechanism is able 
to transport the insolation energy deep enough within a 
planet to significantly affect the planet's evolution (Guil- 
lot & Showman 2002; Burrows et al. 2007). Enhanced 
atmospheric opacity would cause internal heat to be lost 
more slowly, causing a planet's radius to be larger than 
otherwise at a given age (Burrows et al. 2007). In- 
deed, the more highly irradiated planets are thought to 
have enhanced opacity due to species such as gas-phase 
TiO/VO, tholins or polyacetylenes (Burrows et al. 2008; 
Fortney et al. 2008). These upper-atmosphere absorbers 
result in detectable stratospheres (e.g., Knutson et al. 
2009) and prevent incident flux from reaching deep into 
the atmosphere, causing a large day-night temperature 
contrast, which leads to faster cooling (Guillot & Show- 
man 2002). That some planets are not bloated, though 
they are in similar irradiation environments and have 
otherwise similar properties to bloated planets, may be 
due to differences in evolution history or in core mass 
(Guillot et al. 2006; Burrows et al. 2007). 

Currently, the most promising explanation for the large 
radii of some planets is that they were inflated when the 
tidal circularisation of eccentric orbits caused energy to 
be dissipated as heat within the planets (Bodenheimer 
et al. 2001; Gu et al. 2003; Jackson et al. 2008a; Ibgui 
& Burrows 2009). Indeed, Jackson et al. (2008b) found 
that the distribution of the eccentricities of short-orbit 
(a < 0.2 AU) planets could have evolved, via tidal cir- 
cularisation, from a distribution identical to that of the 
farther-out planets. 

The angular momenta of a star and its planets derive 
from that of their parent molecular cloud, so close align- 
ment is expected between the stellar spin and planetary 
orbit axes. When a planet obscures a portion of its par- 
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ent star we observe an apparent spectroscopic redshift 
or blueshift; which we see depends on whether the area 
obscured is approaching or receding relative to the star's 
bulk motion. This manifests as an 'anomalous' radial 
velocity (RV) and is known as the Rossiter-McLaughlin 
(RM) effect (e.g., Queloz et al. 2000a; Gaudi & Winn 
2007). The shape of the RM effect is sensitive to the 
path a planet takes across its parent star, relative to the 
star's spin axis. Thus, spectroscopic observation of a 
transit allows measurement of A, the sky-projected angle 
between the stellar spin and planetary orbit axes. Short- 
orbit, giant planets are thought to have formed just out- 
side the ice boundary and migrated inwards (e.g., Ida 
& Lin 2004). Thus, A is a useful diagnostic for theories 
of planet migration, some of which predict preservation 
of initial spin-orbit alignment and some of which would 
occassionally produce large misalignments. For example, 
migration via tidal interaction of a giant planet with a 
gas disc (Lin et al. 1996; Ward 1997) is expected to 
preserve spin-orbit alignment, whereas migration via a 
combination of planet-planet scattering and tidal circu- 
larisation of a resultant eccentric orbit is able to pro- 
duce a significant misalignment (e.g., Rasio & Ford 1996; 
Chatterjee et al. 2008; Nagasawa et al. 2008). To 
date, A has been determined for 14 systems (Fabrycky 
& Winn 2009; Gillon 2009; Triaud et al. 2009) and for 
3 of those a significant misalignment was found: X0-3b 
(A = 37.3 ±3.7 deg, Winn et al. 2009b; see also: Hebrard 
et al. 2008), HD 80606b (A = 59^11 deg, Gillon 2009; see 
also: Moutou et al. 2009; Pont et al. 2009; Winn et al. 
2009c) and WASP-14b (A = -33.1 ±7.4 deg, Johnson et 
al. 2009; see also: Joshi et al. 2009). 

In this paper, we present the discovery of the transit- 
ing extrasolar planet WASP-17b, which is the least-dense 
planet currently known and the first planet found to be 
in a probable retrograde orbit. 

2. OBSERVATIONS 

WASP-17 is a = 11.6, F6 star in Scorpius. It was 
observed by WASP-South (PoUacco et al. 2006) from 
2006 May 04 to 2006 August 18, again from 2007 March 
05 to 2007 August 19 and again from 2008 March 02 
to 2008 April 19. These observations resulted in 15 509 
usable photometric measurements, spanning two years 
and from two separate fields. A transit search (Collier 
Cameron et al. 2006) found a strong, 3.7-day periodicity 
(Figure [J). 

A full transit of WASP-17 was observed in the /c-band 
with EulerCAM on the 1.2-m Euler-Swiss telescope on 
2008 May 06. The telescope was defocused to give a 
mean stellar profile width of 4". Over a duration of 6 
hours, 181 frames were obtained with a range of expo- 
sure times of 32-98 s — the exposure time was tuned 
to keep the stellar peaks constant. Observations began 
when WASP-17 was at airmass 1.07; the star then passed 
through the meridian before reaching airmass 1.8 when 
observations ended at twilight. The resulting light curve 
and the residuals about the model fits (|j4|) are shown in 
Figure [2] and the photometry is given in Table [1] 

Using the CORALIE spectrograph mounted on the 
Euler-Swiss telescope (Baranne et al. 1996; Queloz et 
al. 2000b), 5 spectra of WASP-17 were obtained in 2007, 
16 more in 2008, and a further 20 in 2009. Three high- 
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Fig. 1. — WASP-South discovery light curve, phase-folded with 
the ephemeris of Table|4] Points with error > 0.05 mag (3 o"„g^j^„) 
were clipped for display. From the WASP discovery photometry we 
found a high probability (0.74) of WASP-17 being a main-sequence 
star and a zero probability of the companion having Rp < 1.5 Rjup 
(Collier Cameron et al. 2007). As such, the system did not ful- 
fill one of our usual selection criteria, 'P{Rp < 1.5Rjup) > 0.2, for 
follow-up spectroscopy (Collier Cameron et al. 2007). We therefore 
advise other transit surveys to exercise caution in rejecting candi- 
dates on the basis of size, so as not to miss interesting systems like 
WASP-17. 
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Fig. 2. — Upper panel: Euler /c-band light curve (red circles) 
taken on 2008 May 06. Overplotted are the best-fitting model 
transits (solid lines; Q from the parameters of Table |4] consult 
that table for the key to the colour and symbol scheme. Lower 
panel: Residuals about the model fits. For the green and the 
magenta models (triangles) the noise is the same (rms = 1140 ppm; 
red noise = 840 ppm — calculated using the method of Gillon et 
al. 2006); for the blue model (squares) the noise is slightly higher 
(rms = 1210 ppm; red noise = 920 ppm, Gillon et al. 2006). The 
mean theoretical error is 800 ppm. 
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TABLE 1 

-BAND PHOTOMETRY OF WASP-17 



H JD-2 450 000 
(days) 


Relative flux 


O'flux 


4592.677426 


0.999803 


0.000959 


4592.678457 


0.998623 


0.000956 


4592.679703 


0.999471 


0.000670 


4592.923564 


1.00053 


0.00141 


4592.926643 


1.00034 


0.00145 



Note. — This table is presented in its en- 
tirety in the electronic edition of the Astrophysi- 
cal Journal. A portion is shown here for guidance 
regarding its form and content. 

precision spectra were obtained in 2008 with the HARPS 
spectrograph (Mayor et al. 2003), based on the 3.6-m 
ESO telescope. RV measurements were computed by 
weighted cross-correlation (Baranne et al. 1996; Pepe 
et al. 2005) with a numerical G2-spectral template. 
RV variations were detected with the same period found 
from the WASP photometry and with semi-amplitude of 
^50 ms^^, consistent with a planetary-mass companion. 
The RV measurements are listed in Table [2] and arc plot- 
ted in Figure [3l 

To test the hypothesis that the RV variations are due 
to spectral line distortions caused by a blended eclipsing 
binary, a line-bisector analysis (Queloz et al. 2001) of 
the CORALIE and HARPS cross-correlation functions 
was performed. The lack of correlation between bisec- 
tor span and radial velocity (Figure |4]), especially for the 
high precision HARPS measurements, supports the iden- 
tification of the transiting body as a planet. 

3. STELLAR PARAMETERS 

The combined CORALIE and HARPS spectra from 
2007-2008, co-added into O.OlA steps, give a S/N of 
~100:1. The stellar parameters and elemental abun- 
dances of WASP-17 were determined using spectrum syn- 
thesis and equivalent-width measurements (Gillon et al. 
2009b; West et al. 2009) and are given in Table El In 
the spectra, the Li i 6708A line is not detected (EW 
< 2mA), giving an upper-limit on the Lithium abun- 
dance of log n(Li/H) -I- 12 < 1.3. However, the effective 
temperature of this star implies it is in the lithium-gap 
(Bohm-Vitense 2004) and so the lithium abundance does 
not provide an age constraint. In determining the pro- 
jected stellar rotation velocity (usini) from the HARPS 
spectra, a value for macroturbulence (umac) of 6 kms~^ 
was adopted (Gray 2008) and an instrumental FWHM of 
0.06 A, determined from the telluric lines around 6300 A, 
was used. A best fitting value of v sini = 9.0 ± 1.5 kms^^ 
was obtained. However, if Wmac is lower than the assumed 
6 kms~^ then wsini would be slightly higher, with a 
value of 11 kms^^ obtained if Wmac is assumed to be 
zero. 

We attempted to measure the rotation period of 
WASP-17 by searching for sinusoidal, rotational mod- 
ulation of the WASP light curve (Hebb et al. 2009), 
as may be induced by a non-axisymmetric distribution 
of starspots. Considering periods of 1.05-30 days, the 
best-fitting period is 24.7 days. However, the amplitudes 
of the phase-folded light curves from each camera from 
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Fig. 3. — Upper panel: Relative radial velocity measurements 
of WASP-17 as measured by CORALIE (red circles). The three, 
colour-coded solid lines are the model orbital solutions (Q based 
on the parameters of Table |4] and incorporate the RM effect. As 
the zero-point offset between HARPS and CORALIE is a free pa- 
rameter in the models, the HARPS measurements are shown once 
per model, with corresponding symbols and colours (Tablc|4]l. The 
centre-of-mass velocities of Table |4] have been subtracted. Lower 
panel: Residuals about the model solutions; consult Table |4] for 
the key to the symbol and colour scheme. 
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Fig. 4. — Bisector span versus relative radial velocity for the 
CORALIE (red circles) and HARPS (blue squares) spectra. Av- 
erages of the centre-of-mass velocities and zero-point offsets of Ta- 
ble |3] were subtracted. Bisector uncertainties equal to twice the 
radial velocity uncertainties have been adopted. The Pearson cor- 
relation coeflicient is 0.19. 
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TABLE 2 

Radial velocity measurements of WASP- 17 
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each season are small (2-8mmag). Assuming spin-orbit 
alignment, with wsini = 9.0kms~^, and using the val- 
ues of stellar radius given in Table [4] (see Q , a stellar 
rotation period of 8.5-11 days is expected. 

We estimated the distance of WASP-17 (400 ± 60 pc) 
using the distance modulus, the TYCHO apparent visual 
magnitude {V = 11.6) and the absolute visual magnitude 
of an F6V star {V = 3.6; Gray 2008); we assumed £^(5- 
V) = 0. 

4. SYSTEM PARAMETERS 

The WASP-South and EulerCAM photometry were 
combined with the CORALIE and HARPS RV mea- 
surements in a simultaneous Markov-chain Monte-Carlo 
(MCMC) analysis (Collier Cameron et al. 2007; PoUacco 



TABLE 3 
Stellar parameters of WASP-17 
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< 1.3 
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11.6 


Distance (pc) 
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et al. 2008). The proposal parameters we use are: Tc, 
P, AF, Ti4, b, Ki, M^, ecosw, esincj, wsinzcosA and 

sin i sin A. Here Tc is the epoch of mid-transit, P is 
the orbital period, AF is the fractional flux-deficit that 
would be observed during transit in the absence of limb- 
darkening, is the total transit duration (from first to 
fourth contact), b is the impact parameter of the planet's 
path across the stellar disc, Ki is the stellar refiex veloc- 
ity semi-amplitude, is the stellar mass, e is the orbital 
eccentricity and uj is the argument of periastron. 

At each step in the MCMC procedure, each proposal 
parameter is perturbed from its previous value by a 
small, random amount. From the proposal parameters, 
model light and RV curves are generated and is cal- 
culated from their comparison with the data. A step is 
accepted if is lower than for the previous step; a step 
with higher x^ may also be accepted, the probability for 
which is lower for larger Ax^- In this way, the param- 
eter space around the optimum solution is thoroughly 
explored. Provided the probability of accepting a step 
of higher x^ is chosen correctly, then the distribution of 
points for an MCMC chain gives the standard errors on 
the parameters (e.g.. Ford 2006). 

We place a prior on that, via a Bayesian penalty, 
causes its values in accepted MCMC steps to approxi- 
mate a Gaussian distribution with mean Mq and stan- 
dard deviation ctm = 0.1 Mq, where Mq is the initial 
estimate of M*. To determine Mq and to estimate the 
star's age, an evolutionary analysis (Hebb et al. 2008) 
was performed. In that, an initial MCMC run was 
used to determine the stellar density, which depends on 
the shape of the transit light curve and the eccentric- 
ity of the orbit. The stellar evolution tracks of Girardi 
et al. (2000) were then interpolated using this stellar 
density and using the stellar temperature and metal- 
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licity from the spectral analysis (Figure O- This sug- 
gests that WASP-17 has evolved off the zero-age main 
sequence, with a mass of 1.20^o }5 -l^o ^^i^ ^S^ of 
S.Otli Gyr. This stellar mass was used as the initial es- 
timate, Mq, in the MCMC solution (Case 1) presented 
in the second column of Table HI The best-fitting ec- 
centricity is non-zero at the 2-a level (e = 0.129lQ'Qgg; 
uj = 290]']}g^ deg). The best-fitting planet radius is large 
but uncertain {Rp = 1.74j;Q23 Rjup). This uncertainty 
results from e and uj being poorly constrained by the RV 
data, causing the velocity of the planet during transit, 
and therefore the distance travelled (i.e. the stellar ra- 
dius), to be uncertain. The planet radius is related to 
the stellar radius by the measured depth of transit, so it 
too is uncertain. 
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Fig. 5. — Modified H-R diagram: inverse cube-root of stellar den- 
sity versus stellar effective temperature. The former quantity is 
purely observational, measured directly from the light curve in ini- 
tial MCMC runs. The latter quantity is derived from the spectral 
analysis. The best-fitting values of the three models for WASP-17 
are depicted using the key of Table |4] Evolutionary mass tracks 
(solid, labelled lines) and isochrones (100 Myr, solid; 1 Gyr, dashed; 
2 Gyr, dotted; 4 Gyr, dot-dashed) from Girardi et al. (2000) for 
[Fe/H] = —0.25 are plotted for comparison. 



The stellar age from the first MCMC solution and 
V sin i from the spectral analysis are consistent with 
WASP-17 being young. As such, a second MCMC anal- 
ysis was performed, with a main sequence (MS) prior on 
the star and with eccentricity a free parameter. With 
the MS prior, a Bayesian penalty ensures that, in ac- 
cepted steps, the values of stellar radius are consistent 
with those of stellar mass for a main-sequence star: the 
probability distribution of i?* has a mean of i?o = Mq '* 
(Tingley & Sackett 2005), where Rq is the initial estimate 
of i?*, and a standard deviation of — 0.8{Ro/Mo)<J]\f. 
An initial MCMC analysis was used to determine stel- 
lar density, which was used as an input to an evolution- 
ary analysis (Figure [5]) . That suggests a stellar mass of 
l.lQtSJ !^^ Mq and a stellar age of 1.2tli Gyr. This stellar 
mass was used as the start value in the MCMC solution 
(Case 2) presented in the third column of Table ID The 
MS prior results in a smaller stellar radius and, there- 
fore, a smaller planetary radius {Rp = 1.51 ± 0.10 Rjup)- 
The MS prior on stellar radius, together with the prior 
on stellar mass, effectively places a prior on stellar den- 
sity, p*, forcing stellar density toward the higher values 



typical of a MS star. Therefore, as 

(1 - e2)3/2 



a more eccentric orbit (e — 0.237tQ;ogg; uj — 

278.0l5 g deg) results. The uncertainties on each of the 
parameters affected by the MS prior are artificially small 
due to the MS prior not taking full account of uncer- 
tainties involved (e.g., in the theoretical mass-radius re- 
lationship). 

As the detection of a non-zero eccentricity in the first 
MCMC solution is of low significance, a third solution 
was generated, with an imposed circular orbit and no MS 
prior. Again, an initial MCMC analysis was performed to 
determine stellar density, which was used as an input to 
an evolutionary analysis (Figure [5]). From that, a stellar 
mass of 1.25 ± 0.08 Mq and a stellar age of S.llg § Gyr 
was found. This stellar mass was used as the start value 
in the MCMC solution (Case 3) presented in the fourth 
column of Tabled] The circular orbit causes the velocity 
of the planet during transit to be higher than in the two 
eccentric solutions. This results in a larger stellar radius 
and, as the depth of transit is fixed by measurement, in 
a larger planet radius {Rp = 1.97 ± 0.10 Rjup). 

For each model, the best-fitting transit light curve is 
shown in Figure [2] and the best-fitting RV curve is shown 
in Figure [H To help decide between the three cases pre- 
sented, a more precise determination of stellar age, stel- 
lar radius or orbtial eccentricity would be useful. It is 
currently difficult to reliably determine the age of stars 
older than 1-2 Gyr (e.g., Sozzetti et al. 2009 and refer- 
ences therein). Stellar radius could be calculated from 
a precise parallax determination. WASP-17's parallax 
is predicted to be 2.5 mas, which will be measurable to 
good precision by the forthcoming Gaia mission (Jordi 
et al. 2006), which is expected to achieve an accuracy of 
7 /ias a.t V — 10. Eccentricity can be better determined 
using a combination of two methods: (i) Take a number 
of high-precision RV measurements (which best constrain 
e sinw), focusing on those phases at which the differences 
between the models are greatest (Figure [3]). (ii) Observe 
the secondary eclipse; the time of mid-eclipse constrains 
e cos Lu and the eclipse duration more weakly constrains 
esinw (Charbonneau et al. 2005). 

We adopt Case 1 as our preferred solution; we note 
that if WASP-17 proves to be young then Case 2 will 
be indicated, and Case 3 will be indicated if the planet's 
orbit is found to be (near-)circular. 

4.1. A retrograde orbit? 

The RM effect was incorporated in the MCMC anal- 
yses with free parameters v sin i cos A and v sin i sin A 
(Figure [H Table |4]). The three RV measurements 
during transit suggest a large spin-orbit misalignment 
(A ~ — 150deg), indicating that the planet is orbiting 
in a sense counter to that of stellar rotation. The RV 
RMS about the fitted model is 31.6ms~^. Comparing 
a borderline prograde-retrograde orbit (A = — 90deg, 
wsini = 9.0kms~^, b = 0.355), the first in-transit point 
is discrepant by 5.0 cr and the RV RMS is 40.2 ms~^ 
(Figure [6]). Comparing a prograde orbit (A = Odeg, 
wsini = 9.0kms~^, b = 0.355), the first and third in- 
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transit points are discrepant by 6.0 a and 3.8(7 respec- 
tively, and the RV RMS is 45.7 ms'^ (Figure [6]). 

The fitted amplitude of the RM effect suggests v sin i « 
20kms~^, which is higher than determined in the spec- 
tral analysis (Table [3]). This could be because the am- 
plitude of the RM effect is currently liable to be over- 
estimated (Winn et al. 2005; Triaud et al. 2009), due 
to the manner in which the RVs are extracted from the 
spectra. At present, the effective velocity of the spectro- 
scopic cross-correlation function (CCF) is measured by 
fitting a gaussian. However, at values of wsini signifi- 
cantly greater than the intrinsic width of the CCF for 
a slowly-rotating star, the travelling bump in the profile 
that is the spectral signature of the planet's silhouette 
becomes partially resolved (Gaudi & Winn 2007). The 
CCF profile will become slightly asymmetric when the 
planet is near the limb, and this may bias the velocity 
measured by gaussian-fitting to a greater value than the 
RV of the centroid of the unobscured parts of the star. 




-150 ^ — ■ — ■ — ' — ' — ■ — ■ — ■ — ' — ■ — ■ — ■ — ' — ^ 
0.95 1 1.05 

Orbital phase 



Fig. 6. — Zoom-in on the spectroscopic transit region of Figure |3] 
showing the fitted RM effect more clearly. The red circles are 
CORALIE measurements and the coloured lines are the best-fitting 
models (Table [Jjl . We show for comparison the RM effects that 
would result from perpendicular (A = — 90deg; black, dot-dashed 
line) and aligned (A = Odeg; black, dashed line) spin-orbit axes; in 
both cases vsini = 9.0kms~^ and b = 0.355 were fixed. 



4.2. Transit times 

We measured the times of WASP-17b's transits to 
search for transit timing variations, as may be induced 
by a third body (e.g., Holman & Murray 2005; Agol et al. 
2005). The model hght curves were stepped in time over 
the phometric data around the predicted times of transit, 
and was calculated at each step. The times of mid- 
transit were found by measuring the minima and the 
uncertainties were determined via bootstrapping. The 
calculated times of mid-transit, Tc, and the differences, 
O — C, between those times and the predicted times, as- 
suming a fixed epoch and period (Table 2]), are given in 
Table [S] No significant departure from a fixed ephemeris 
is seen. 



5. DISCUSSION 

WASP- 17b is the least dense planet known, with a den- 
sity of 0.06-0.14pj. TrES-4, the previous least dense 
planet, has a density of 0.15pj (Sozzetti et al. 2009), 
and HD 209458b, the most studied transiting planet, has 
a density of 0.27/9j (Torres et al. 2008). WASP-17b's 
radius of 1.5-2 Rj up is larger than predicted by standard 
planet evolution models. For example those of Fortney et 
al. (2007) imply a radius of at most l.SRjup (the value 
for a 1-Gyr-old, coreless planet of 0.41 Mjup, receiving 
more stellar flux, at a distance of 0.02 AU, where each of 
these values errs on the side of inflating the radius). 

Burrows et al. (2007) showed that an enhanced atmo- 
spheric opacity can delay radius shrinkage, leading to a 
larger-than-otherwise planet radius at a given age. En- 
hanced opacities may result from super-solar metallicity, 
the presence of clouds/hazes, or the effects of photol- 
ysis or non-equilibrium chemisty. One might expect a 
low planetary atmospheric opacity due to the sub-solar 
metallicity of the WASP-17 star: [Fe/H] = -0.25. How- 
ever, as WASP-17b is highly irradiated, its atmospheric 
opacity is expected to be high due, for example, to the 
presence of TiO and VO gases (Burrows et al. 2008; 
Fortney 2008). Ibgui & Burrows (2009) found that an 
atmospheric opacity of 3 x solar is sufficient to account 
for the radius of HD 209458b. The very large radius of 
WASP-17b and its moderate age suggest that enhanced 
opacity alone, even of 10 x solar, is insufficient to account 
for the planet's bloatedness. 

It has been proposed (Bodenheimer et al. 2001; Gu 
et al. 2003; Jackson et al. 2008a) that tidal dissipa- 
tion associated with the circularisation of an eccentric 
orbit is able to substantially inflate the radius of a short- 
orbit, giant planet. If a planet is in a close (a < 0.2 AU), 
highly eccentric (e > 0.2) orbit then planetary tidal dissi- 
pation will be significant and will shorten and circularise 
the orbit. Orbital energy is deposited within the planet 
interior, leading to an inflated planet radius. This pro- 
cess is accelerated by higher atmospheric opacities: as 
the planet better retains heat, shrinking of the radius 
is retarded, and a larger radius causes greater tidal dis- 
sipation. Higher eccentricities result in stronger tides, 
and thus in greater tidal dissipation. The rate at which 
energy is tidally dissipated within a body is inversely pro- 
portional to its tidal quality factor, Q' , which is the ratio 
of the energy in the tide to the tidal energy dissipated 
within the body per orbit (e.g., Ogilvie & Lin 2007). 

Ibgui & Burrows (2009) created a tidal dissipation 
model and applied it to HD 209458b, which is bloated 
to a lesser degree than WASP-17b. A custom fit is re- 
quired to find possible evolution histories for the WASP- 
17 system, but the similarity of HD 209458 (compare 
Tables |3] and [3] from this paper with Table 1 in Ibgui 
& Burrows (2009) and references therein) permits com- 
parison. Ibgui & Burrows' (2009) HD 209458b model 
suggests that tidal heating could produce even WASP- 
17b's maximum likely radius (i?p w 2Rjiip) if, for exam- 
ple, it evolved from a highly eccentric (e « 0.79), close 
(a « 0.085 AU) orbit, with moderate tidal dissipation 
(Qp « 10^-^^; Q'^ w 10''°) and solar atmospheric opac- 
ity. The final semimajor axis of this particular model 
was shorter than that of WASP-17b, but within 10 per 
cent. Such an eccentric, short orbit seems reasonable as 
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TABLE 4 

System parameters of WASP-17 



Parameter 



Case 1 (adopted) 



Case 2 



Case 3 



Extra constraints 
Graph colour 
Graph symbol 
Graph line 
Stellar age (Gyr) 



blue 
squares 
solid 

O Q + 0.9 



MS prior 
green 
upwards triangles 
dashed 



1.2 



+2.8 



e = 
magenta 
downwards triangles 
dotted 

n -,+1.1 
■^•-^-0.8 



P(d) 

Tc (HJD) 
Tl4 (d) 
Ti2 = T34 (d) 
AF = H2 /r2 
b = a cos i I Rt, 



o 70^4/11 7+0.0000072 
o. lOdl^il I _o. 0000073 



2454559.18102 



+0.00028 
0.00028 



0.1822 
0.0235 



+0.0019 
-0.0023 
+0.0019 
-0.0030 



01672+° °°°^^ 



0.352 



+0.075 
0.316 



S 7QK441 7+0.0000073 
o. I Odl^il / -0.0000074 



2454559.18100 



+0.00027 
0.00028 



n I09c:+0.0017 
U-l»^O_0.0017 

U.UZJD_Q QQJ^g 

01674+°-°°°2'' 



0.355 



+0.068 
0.111 



3.7354414 
2454559.18096 



+0.0000074 
0.0000074 
+0.00028 
0.00028 



0.1824 



+0.0016 
-0.0016 



n23q+°°°" 

0.01b78_Q 150027 
„ ,~„+0.084 

U.O(U_o Qgg 



Ki (kms-1) 
7 (kms"-') 

7HARPS -7COR. (kms-l) 
a (AU) 



n n56q+° °°^^ 

U.UODy_o QQ53 
-49 5128+° °°!'^ 
0267+°-°°3'' 
0501+° °°" 

U.UOUi_(J QQJg 



n5q2+° °°5** 
-49 5128+°-°°^* 
0.0297t°:°°f3 
0.0494+°°°" 



0564+° °°'5^ 

U.U0D4_Q 005J 

-49 5125+° °°" 
U.UZOCi_o oojg 
0507+°-°°" 



i (deg) 
e cos Lu 
e sinuj 
e 

uj (deg) 

■^mid-cclipsc 

T58 (d) " 

756 = Tyg (d) >> 



87.S 



+2.0 
1.0 



88.16 



036+° °^" 



-0.10 
0.129 
290 



+0.13 
0.13 
+0.106 
0.068 
106 
16 



n t;9q+0.021 
u.jzo_o 020 

152+°°'*° 
U.10Z_o 034 

"•0100-0.0046 



+0.58 
0.45 



86.95 



034+° °25 
u.uoi_Q Q24 

r, r,oo+0.071 
U.ZOO_o Q7Q 

+0.068 
0.069 
+8.2 
5.6 

n ^99+0.016 

(, ,,-+0.017 
'J-J^J^' -0.015 

0140+° °°22 



0.237; 
278.0 



+0.87 
0.63 



0. (fixed) 



0.5 (fixed) 
-0.0016 (fixed) 
0.0239to:ol5l7 (fixed) 



0.1824+°°°« 



A (deg) 

vsini (kms"^) 



-147; 
20.0 



+49 
11 
+69.2 
5.2 



-148.7; 



19.1 



9.3 
+5.7 



-149.3: 



19.1 



+5.3 



(Mo) 
R* (R©) 

logs* (cgs) 



1.20 
1.38 
4.23 
0.45 



+0.12 
0.12 
+0.20 

0.18 
+0.12 
0.12 
+0.23 
0.15 



1 16+°-l2 
^■^"-0.12 

1 200+°°*! 

l.ZUU_Q QgQ 

4 341+° °68 
^•"-'^^-0.068 

67+°-l'5 

-0.13 



, 9C-+0.13 
1-^3-0.13 

l-ODD-0.073 
4 14Q+O.O32 
^■^^•J-0.031 

n 090+0.035 

U.JZO_g Q2g 



Mp (Mjup) 
Rp (R,jup) 
loggp (cgs) 

PP (Pj) 
7p,A=o (K) 



490+° °53 
u.tyu_Q Q5g 



1.74 
2.56 



+0.26 
0.23 
+0.14 
-0.13 

092+° °^* 

u.uyz_Q Q32 
1662l}}3 



0.496 
1.51 



+0.064 
0.060 
+0.10 
0.10 



9 006+0-086 
z.uyo_Q Qg3 

U.1^4-0.031 

155711 



498+0.059 
u.'iye5_Q Qjg 



+0.10 
0.10 
+0.051 



1.97 
2.466 



n nfi/i8+° °i°6 
0.0b4S_(j oijgp 

1756+^^ 



Three solutions are presented (Cases 1, 2 and 3), each with different constraints as described in the text (Ejl. 
" T58: total eclipse duration. Tsg = T78: eclipse ingress/egress duration. 
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TABLE 5 
Transit times 









ij — o 




/"H" T"^»^ 

\n..}D ) 


(min) 


(iQin) 


— 179 


245o»yu.54»5U 


6.2 


16.6 


— 175 


245ci9Uo.48l9o 


5.5 


4.7 


— 171 


z45o9zU.4zz4o 


3.6 


2.8 


— -Luy 




O.U 


— 12 2 


-96 


2454200.57081 


4.4 


-11.2 


-92 


2454215.52194 


2.7 


2.3 


-77 


2454271.55701 


4.1 


7.2 


-73 


2454286.49409 


8.3 


0.5 


-69 


2454301.45133 


8.2 


22.8 


-61 


2454331.32310 


9.3 


5.8 


-1 


2454555.43662 


6.4 


-12.9 


2 


2454566.65046 


8.3 


-2.0 


9 


2454592.80049 


0.55 


0.79 



planets in highly eccentric, quite short orbits are known: 
HD 17156b (e = 0.676, a = 0.162 AU; Winn et al. 
2009d), HD 37605b (e = 0.74, a = 0.26 AU; Cochran et 
al. 2004), HD 80606b (e = 0.934, a = 0.45 AU; Moutou 
et al. 2009). As for the strength of tidal dissipation, 
Jackson et al. (2008b) found similar best-fitting values 
(Q^ = 10^-^; Q'^, = 10^-^) when matching the current 
eccentricities of short-orbit (a < 0.2 AU) planets with 
the eccentricities of farther out planets, from which they 
presumably evolved. 

The limited radial-velocity measurements during tran- 
sit give a strong indication that WASP- 17b is in a retro- 
grade orbit. As the angular momenta of a star, its proto- 
planetary disc, and hence its planets, all derive from that 
of the parent molecular cloud, WASP- 17b presumably 
originated in a prograde orbit. As a gas giant, WASP-17b 
is expected to have formed just outside the ice boundary 
(~3 AU) and migrated inwards to its current separation 
of 0.05 AU (e.g., Ida & Lin 2004). Migration of a giant 
planet via tidal interaction with a gas disc is expected 
to preserve spin-orbit alignment (Lin et al. 1996; Ward 
1997) and is thus unable to produce a retrograde orbit. 
Alternatively, migration via a combination of star-planet 
scattering (Takeda et al. 2008) or planet-planet scat- 
tering (e.g., Rasio & Ford 1996; Chatterjee et al. 2008; 
Nagasawa et al. 2008) and tidal circularisation of the 
resultant eccentric orbit is able to produce significant 
misalignment. 

In addition to inclined orbits, scattering is able to pro- 
duce highly eccentric orbits (e.g.. Ford & Rasio 2008), 
which have been found to be common and are neces- 
sary if tidal circularisation is to inflate planetary radii 
(e.g., Jackson et al. 2008b), whereas planet-disc interac- 
tions seem unable to pump eccentricities to large values 
(e > 0.3; D'Angelo et al. 2006). Nagasawa et al. (2008) 
carried out orbital integrations of three-planet systems: 
three Jupiter-mass planets were initially placed beyond 
the ice boundary (5AU, 7.25 AU, 9.5 AU) in circular 
orbits, with small inclinations (0.5°, 1°, 1.5°), around 
a solar-mass star. They found that a combination of 
planet-planet scattering, the Kozai mechanism (the os- 
cillation of the eccenticity and inclination of a planet's or- 
bit via the secular perturbation from outer bodies; Kozai 



1962) and tidal circularisation produces short-orbit, gi- 
ant planets in ~30% of cases. The Kozai mechanism 
is most effective when the scattered, inner planet has 
an inclined orbit. A broad spread was seen in the in- 
clination distribution of the short-orbit planets formed, 
including planets in retrograde orbits. Therefore, we sug- 
gest that WASP- 17b supports the hypothesis that some 
short-orbit, giant planets are produced by a combina- 
tion of scattering, the Kozai mechanism and tidal circu- 
larisation. The observation of the RM effect for more 
short-orbit planets is required to measure the size of the 
contribution. 

For planet-planet or star-planet scattering to have in- 
fluenced WASP-17b's orbit in the past, one or more stel- 
lar or planetary companions must have been present in 
the system. Sensitive imaging could probe for a stel- 
lar companion and further radial velocity measurements 
are necessary to search for stellar or planetary compan- 
ions. It will be worthwhile looking for long-term trends 
in the RVs to detect farther-out planets that might have 
been involved in past scattering. A straight-line flt to the 
residuals of the RV data about the model flts indicates 
no significant drift over a span of 622 days (e.g., for Case 
1 the drift is — 17± 11 ms^^). In their three-planet in- 
tegrations, Nagasawa et al. (2008) found that in 75% of 
cases one planet is ejected, a planet collides with the host 
star in 22% of cases, and two planets are ejected in 5% 
of cases. They also found that, since a small difference 
in orbital energy causes a large difference in semimajor 
axis in the outer region, the final semimajor axes of outer 
planets are widely distributed (peak at ~15AU, with a 
large spread). Therefore, it is possible that WASP-17 
is now the only giant planet in the system or that the 
outer planets are in long orbits, which are difficult to 
detect with the RV technique. 

The discovery of WASP-17b extends the mass-radius 
distribution of the 62 known transiting exoplanets (Fig- 
ure [7]). WASP- 17b has the largest atmospheric scale 
height (1100-2100 km) of any known planet by up to 
a factor 2, due to its very low surface gravity and mod- 
erately high equilibrium temperature. The ratio of pro- 
jected areas of planetary atmosphere to stellar disc of 
WASP-17b is 1.9-2.7 times that of HD 209458b and 2.4- 
3.4 times that of HD 189733b, for both of which success- 
ful attempts at measuring atmospheric signatures have 
been made (e.g., Charbonneau et al. 2002, Desert et al. 
2009). Thus, although WASP-17 is fainter and has a 
larger stellar radius, the system is a good prospect for 
transmission spectroscopy. 
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